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Abstract 

We consider a spherical thick shell immersed in two different spherically symmetric 
space-times. Using the fact that the boundaries of the thick shell with two embed- 
ding space-times must be nonsingular hypersurfaces, we develop a scheme to obtain 
the underlying equation of motion for the thick shell in general. As a simple exam- 
ple, the equation of motion of a spherical dustlike shell in vacuum is obtained. To 
compare our formalism with the thin shell one, the dynamical equation of motion 
of the thick shell is then expanded to the first order of its thickness. It is easily 
seen that the thin shell limit of our dynamical equation is exactly that given in the 
literature for the dynamics of a thin shell. It turns out that the effect of thickness 
is to speed up the collapse of the shell. 
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I. Introduction 



The thin shell formalism of general relativity has found wide applications in general 
relativity and cosmology [1- 3] . Studies on gravitational collapse, dynamics of bubbles 
and domain walls in inflationary models, wormholes, signature changes, structure 
and dynamics of voids in the large scale structure of the universe are some of the 
applications. Thin shells are considered as idealized zero thickness objects, with 
a 5-function singularity in their energy-momentum and Einstein tensors. This is 
regarded to be an idealization of a real shell with a finite thickness. However, the 
dynamics of a real thick shell has been rarely discussed in the literature because of 
the complexity one is faced with when trying to define it within general relativity 
and to find its exact underlying dynamical equations. The outstanding paper that 
modifies the Israel thin shell equations to treat the motion of spherical and planar 
thick domain walls is that of Garfinkle and Gregory [4]. Their work deals with an 
expansion of the coupled Einstein-scalar equations in powers of the thickness of the 
domain wall (see also [5]). According to the results of that paper, the effect of 
thickness in the first approximation is to reduce effectively the energy density of the 
wall compared to the corresponding thin domain wall, and therefore to increase the 
collapse velocity of the wall in vacuum. 

In this paper we first generally suggest the proper matching conditions on the bound- 
aries of the spherical thick shell embedded in an inner and an outer spherically sym- 
metric space-time. As a simple example, these matching conditions are then used 
to investigate the motion of spherical dustlike thick shell in vacuum. 

II. The junction conditions 

Consider a spherically symmetric thick shell with two boundaries Si and S2 dividing 
the space-time into three-regions: for inside the inner boundary Si,A4oMt for 

outside the outer boundary E2, and M. for the thick shell having two boundaries 
El and E2. First of all, let us write down the appropriate junction condition on 
each boundary E^ (j=l,2) treated as a (2-|-l)-dimensional timelike hypersurface. 
We expect the continuity of the second fundamental form of Ej, or the extrinsic 
curvature tensor Kab of Ej, so that we can consider Ei(E2) as a boundary surface 
separating Ai region from A4in {-Mout) ■ This crucial requirement is formulated as 



where the square bracket indicates the jump of Kab across E^, Latin indices range 
over the intrinsic coordinates of E^ denoted by {Tj,9, (p), where Tj is the proper time 
of Ej. In particular, the angular component of Eq. (1) on each boundary is written 
as 
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where the superscript +(— ) refers to the side of towards which the corresponding 
unit spacehke normal vector n"(— n") points. This means that on Ei(S2), the 
superscript + refers to the region M.{M.out) and the superscript — refers to the 
region Aiin{M.). Adding Eqs. (2) and (3), we obtain equation 
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In the next section we will apply the general equation (4) to the special case of a 
collapsing spherical dust shell in vacuum. 



III. Collapse of a spherical dust thick shell in vacuum 

Consider the Lemaitre-Tolman-Bondi (LTB) metric to describe the dust thick shell. 
In the synchronized comoving coordinates (r, r, 9, (p) the metric is written in the 
form [6] 

ds^ ^-dT'' + Y^^^dr^ + {r, t) {d9^ + sin^ Odip^) , (5) 

where overdot and prime denote partial differentiation with respect to r and r, 
respectively, and E{r) is an arbitrary real function such that E{r) > —1. Then the 
corresponding Einstein field equations turn out to be 

R\r,T)^E{r) + ^, (6) 
F'(r) 

where p{r,T) is the energy density of the matter fluid in M., and F{r) is another 
arbitrary real smooth function such that F{r) > 0. Furthermore, in order to avoid 
shell crossing of dust matter during their radial motion, we require R'{r, r) > 0. 
This together with the assumption of positive mass density p{r, r) > 0, implies that 
F'(r) > 0. The induced intrinsic metric on may be represented as 



ds'^ 



= -drf + R^ATj){de^ + sin^ edcp^) {j = 1, 2), (8) 



where Rjijj) being the proper radius of S^. For simplicity, we may assume that 
the boundaries Si and E2 are comoving with respect to the LTB geometry. This 
requires that the peculiar velocity of E^ measured by the comoving observers of Ai 
to be zero. Then the matching relations yield 

T = Tj + constant {i = 1,2). (9) 

Now, we may define the constant comoving thickness of the shell as follows 

25 = r2-ri, (10) 
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where ri and r2 are comoving radii of the boundaries Ei and S2, respectively. We 
assume now the spherical dust thick shell to be immersed in vacuum. In this case, 
the space-time exterior to the shell is Schwarzschild, and the interior is taken to be 
Minkowski flat space-time. Different terms appearing in the equation (4) may now 
be explicitly derived. Using the metric (5) together with Eqs. (6), (7) and (9), we 
compute the relevant extrinsic curvature tensors in the region Ai as 
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where Rj = R{rj,T). Furthermore, the following expressions for the relevant extrin- 
sic curvature tensors in Aiin and A4out also hold [3] 
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where TZ{r2) is the Schwarzschild radius of the spherical shell within the comoving 
surface r2. Now, to obtain the dynamical equation of the thick shell, we first expand 
the following quantities in a Taylor series around tq, the mean comoving radius of 
the thick shell: 



R{rj,T) = i?(ro, r) + e,-5i?'(ro, r) + 0{6''), 



(13) 



F(rj) = F(Ro) + ej5F'{ro) + 0(5''), 



(14) 
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where ei = —1 and 62 = +1. Using Eqs. (13), (14) and (15) in the expressions (11) 
and (12) and keeping only terms up to the first order of 5, we obtain 
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where Rq = R{ro,r). Substituting Eq. (16) into Eq. (4) and noting that for the 
metric LTB, F(ro) plays just the role of the Schwarzschild radius of the part of the 
spherical thick shell within the comoving surface ro denoted by 7?.(ro), we obtain 
after some rearrangement the dust thick shell's equation of motion written up to 
the first-order in S: 



F'(ro) 
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with 



, - « + RR'.i^) + ^(T^Vo) + . (17) 

a = \/l + fi§ , /J=^1 + BS-^. (18) 

This is the generalization of thin shell dynamical equation up to the first order of 
the thickness. It is now interesting to verify the thin shell limit of this thick shell 
dynamical equation. In order to do this, consider the following definition for the 
surface energy density of the infinitely thin shell[l]: 

p{r,T)dn, (19) 

where n is the proper distance in the direction of the normal and 2e is the physical 
thickness of the shell. With the metric (5), Eq. (19) takes the form 

/ P(r,r) ^ ' [ dr. (20) 
J-5 Jl + Eir) 



E{r) 

Using Eqs. (6) and (7), we find that Eq. (20) can be written as 

SnGa^ / ^ ^ -.dr- (21) 

J-^ /?yi + i?2_£M 

Using Eqs. (13) and (14), we may integrate Eq. (21) up to the first order in 5 to 

get 

S7rGa^2S ^£^k^=== + O(d^). (22) 

Substituting (22) into Eq. (17), we get the following result 
a- P ^ AnGaRo 
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Note that in the zero thickness hmit of the shell, as 6 — > 0, the second term on the 
right hand side of Eq. (23) is regular and goes to zero such that Eq. (23) reduces to 
the Israel's equation of motion for the dust thin shell in vacuum, as given by Israel 
[1], Ipser-Sikivie [7] , Sato[8], and Khorrami-Mansouri [9]. 

To see the explicit effect of the thickness on the dynamics of the thick shell we 
rewrite Eq. (23) as 

a -13 = AirGaRo, (24) 
where a, the effective surface density, is defined by 

* = -sk (|(" - « ^ *^'°(^) + ^C^'*^"* ^ f '^<'°») ■ <^^' 

We see that Eq. (24) has the same form as the well-known Israel's equation for a 
thin shell with the effective surface density a. Now let us have a closer look at the 
terms within the brackets. Note that for a dust shell starting its collapse at rest, 
the velocity Rq is negative during the collapse, it also becomes more negative with r 
so that R'o < 0, so the combination of RoR'o must be positive. On the other hand, 
the Schwarzschild radius of the shell layers is increased with r so that TZ'{ro) > 0. 
Therefore all the terms within the bracket on the right hand side of Eq. (25) are 
positive. This leads to the result a < a. Now, solving Eq. (24) for we find 

From Eq. (26), it follows that for a given Rq and TZ{ro) as long as Rq > 7^(ro), 
smaller a{d- < a) leads to a larger R^. Therefore the first-order thickness corrections 
to the Israel thin shell approximation speed up the collapse of dust shell in vacuum. 



IV. Conclusion 

We have presented a simple procedure to investigate the dynamics of a spherical 
thick shell embedded in an otherwise spherically symmetric space-time, based on the 
Darmois conditions satisfying on the shell boundaries. As the simplest nontrivial 
example, we applied our scheme to the case of the collapse of a thick shell of dust 
in vacuum and obtained the zero thickness limit of our formalism which is just the 
Israel thin shell equation. It has been shown that the effect of thickness, up to the 
first order in the shell thickness, is to speed up its collapse in vacuum. 
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